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Abstract 

By considering delta function sextupole, octupole, and decapole perturbations and using 
difference action-angle variable equations, we find some useful analytical formulae for the 
estimation of the dynamic apertures of circular accelerators due to single sextupole, single 
octupole, single decapole (single 2m pole in general). Their combined effects are derived 
based on the Chirikov criterion of the onset of stochastic motions. Comparisons with nu- 
merical simulations are made, and the agreement is quite satisfactory. These formulae have 
been applied to determine the beam-beam limited dynamic aperture in a circular collider. 

1 Introduction 

One of the preoccupations of the circular accelerator designer is to estimate the influence of 
nonlinear forces on the single particle's motion. These nonlinear forces manifest themselves as 
the systematic and random errors from optical elements, the voluntarily introduced functional 
ones, such as the sextupoles used for the chromaticity corrections and the octupoles used for 
stabilizing the particles' collective motion, or from nonlinear beam-beam interaction forces. 
Even though the nonlinear forces mentioned above compared with the linear forces are usually 
very small, what is observed in reality, however, is that when the amplitudes of the transverse 
oscillation of a particle are large enough, the transverse motion might become unstable and 
the particle itself will finally be lost on the vacuum chamber. Apparently, the above implied 
maximum oscillation amplitudes, A xy , corresponding to the stable motions are functions of the 
specific longitudinal position, s, along the machine, and these functions A Xt y(s) are the so-called 
dynamic apertures of the machine. A reasonably designed machine should satisfy the condition 
A x ,y(s) > M xy (s), where M Xty (s) are the mechanical cross-section dimensions of the vacuum 
chamber. 

Needless to say the dynamic aperture problem in circular accelerators is one of the most 
challenging research topics for accelerator physicists, and the relevant methods adopted to treat 
this problem are quite different analytically and numerically. In this paper we will show how 
single sextupole, single octupole, and single decapole (single 2m pole in general) in the machine 



limit the dynamic aperture and what is their combined effect if there are more than one non- 
linear element. From the established analytical formulae for the dynamic aperture one gets the 
scaling laws which relate the nonlinear perturbation strengths, beta functions, and the dynamic 
apertures. To test the validity of these formulae, comparisons with some numerical simulation 
have been made. As an interesting application the beam-beam limited dynamic apertures in a 
circular collider have been discussed. 



2 Hamiltonian Formalism 

The Hamiltonian formulation of dynamics is best known not only for its deep physical and 
philosophical inspiration to the physicists, but also for its technical convenience in solving various 
nonlinear dynamical problems. The general Hamiltonian for a particle of rest mass, mo, and 
charge, e, in a magnetic vector potential, A, and electric potential, is expressed as: 

2W2 



H(q,p,t) = e^ + c^{p-eAy + m c z j (1) 

where c is the velocity of light and p is the momentum with its components, pi, conjugate to the 
space coordinates, The equations of motion can be readily written in terms of Hamilton's 
equations: 

dpi = &H 

dt d qi { ' 

dt dpi 

For our specific dynamic problems in a circular accelerator it is convenient to chose curvilinear 
coordinates instead of Cartesian ones for us to describe the trajectory of a particle near an 
a priori known closed orbit. The Hamiltonian in the new system, (x, s, y) (where x, s, and y 
denote the coordinates in the Frenet-Serret normal, tangent, and binormal triorthogonal and 
right-handed coordinate system) , is given by Q ~~ Jjfl : 



H t (q,p, t) = e$ + c [ {p x - eA x f + (p y - eA y ) 2 + ( P ° - *^° ) + m c 2 ) (4) 



where p is the radius of curvature and the torsion of the closed orbit is everywhere zero. Since 
it is useful to use the variable, s, as the independent variable rather than time, t, one gets the 
new Hamiltonian by using a simple canonical transformation: 

/ 1 \ 1/2 

H s = -eA s -{l + x/p) l-(E 2 -mlc A )-{p x - eA x f - {p y - eA y f\ - e$ (5) 

Noting that the term \{E 2 — mgC 4 ) in eq. [| is equal to P 2 with P being the total mechanical 
momentum of the particle, by using another trivial canonical transformation: 



- - - Px.y u H s 

q = Q,s = s,p x , y = ~^-,H = — (6) 
Pq Pq 



one gets another Hamiltonian: 



H = - (1 + x/p) (—-(p x - —f - (p v - e -^L) 2 X' 2 - — (7) 

Po K ' ' \Po y Pa } K v Po ' J Po ^ J 
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where Po is the mechanical momentum of the reference particle, and P = Pq + AP. Inserting 
eA s /P in eq. by: 

one gets finally the Hamiltonian which serves as the starting point of most of the dynamical 
problems concerning circular accelerators: 

X 2 B y \ X = 0j y = Q 1 ^ 1 d'^By ( ~ J_ j „ ,\Tl 

2,p Bq Bop ^ nl ox 11 1 

-^ i+ £>-£)'-(*-£)'f-* (9) 

where Bq is the bending magnetic field on the orbit of the reference particle, and B y in general 
is a complex variable. 



3 Analytical formulae for dynamic apertures 

To start with we consider the linear horizontal motion of the reference particle (no energy 
deviation) in the horizontal plane (y=0) assuming that the magnetic field is only transverse 
{A x = Ay = 0) and has no skew fields, and $ is a constant. The Hamiltonian can be simplified 

as 

do) 

where x denotes normal plane coordinate, p = dx/ds, and K(s) is a periodic function satisfying 
the relation 

K(s) = K(s + L) (11) 
where L is the circumference of the ring. The solution of the deviation, x, is found to be 

x = \fe x /3 x (s) cos(0(s) + 4> ) (12) 

where 

(13) 



10 Px(s) 

As an essential step towards further discussion on the motions under nonlinear perturbation 
forces, we introduce action-angle variables and the Hamiltonian expressed in these new variables: 

f s ds' 

•i5M + * (14> 



2 2/3.. 



J= * = =w)[* + (15> 



ff «*) = SM (16> 

Since H ( J, = J/(3 x (s) is still a function of the independent variable, s, we will make another 
canonical transformation to freeze the new Hamiltonian: 

_ 2ttu f s ds' 
L Jo Px(s') 
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Jl = J 

2-kv 



L 



-Ji 



(18) 
(19) 



Before going on further, let's remember the relation between the last action-angle variables and 
the particle deviation x: 



2 J,. },.(.«) cos ( * a - + 



ds' 
o M^) 



(20) 



Being well prepared, we start our journey to find out the limitations of the nonlinear forces on 
the stability of the particle's motion. To facilitate the analytical treatment of this complicated 
problem we consider at this stage only sextupoles and octupoles (no skew terms) and assume 
that the contributions from the sextupoles and octupoles in a ring can be made equivalent to a 
point sextupole and a point octupole. The perturbed one dimensional Hamiltonian can thus be 
expressed: 



TT p z K(s) 2 

H = 1 — T 4- 

2^2^ 3\Bp dx 2 



1 d 2 B z 3t ~ , r . 1 d 3 B z , 



5 L ${s-kL) + 



k=—oc 



4\Bp dx 3 



L ]T 5{s-kL) 



k=—oo 



Representing eq. 21 by action-angle variables (J\ and \&i), and using 



one has 



B z = B (l + xbi + x 2 b 2 + x 3 b 3 
{2J l f3 x (s 1 )f/ 2 



(21) 



(22) 



L 3p 

W X (S2)) 2 



b 2 L cos 3 ^ <>{s ~ kL) 

k=—oo 



+ 



b%L cos 



(23) 



r k=—oo 

where s\ and s 2 are just used to differentiate the locations of the sextupole and the octupole 
perturbations. By virtue of Hamiltonian one gets the differential equations for and J\ 



dJi 
ds 



dHi 



0*i 

ds dJ\ 
d,h (2J l (3 x (s 1 )f/ 2 , dcos 3 ^! 



(24) 
(25) 



ds 



3p 



-b 2 L- 



<Nfi 



E 8{s-kL) 



k=—oo 



± Hs-kL) 



p 



(26) 



fc=— oo 



d*! 2nv , V2J^(3 x ( Si r/ 2 u T 3 ~ 

= 1 i b 2 L cos Wi > d(s — kL) 

ds L p , ' 

' fc=— oo 

2/?2( S2 ) 



+ 



J163 L COS 

fc=— oo 



(27) 



We now change these differential equations to the difference equations which are suitable to 
analyse the possibilities of the onset of stochasticity HQ. Since the perturbations have a 
natural periodicity of L we will sample the dynamic quantities at a sequence of Sj with constant 
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interval L assuming that the characteristic time between two consecutive adiabatic invariance 
breakdown intervals is shorter than L/c. The differential equations in eqs. 2£ and 27 are reduced 
to 

T x =7T(*i, Ji) (28) 

W~ 1 =W 1 (^ 1 ,J 1 ) (29) 

where the bar stands for the next sampled value after the corresponding unbarred previous 
value. 

— (2J!/3 x ( Sl )) 3 / 2 , dcOS 3 ^l (JMS 2 ))\ r dcOS 4 *l 
J\ = J\ 02 L — 03-^" 



3p 



) 3/2 Ji 



1/2 



-b 2 L cos 6 *i + 



2^(s 2 )V 



Ji6 3 Lcos 4 ^ 



(30) 



(31) 



Eqs. 3C and 31 are the basic difference equations to study the nonlinear resonance and the onset 
of stochasticities considering sextupole and octupole perturbations. By using trigonometric 
relation 

m 1 

ETYl 
— — cos(n - m + 2r)9 (32) 
n (m-r)!r! 



cos m 9 cos n9 



one has 



cos 4 (9 



cos 3 9 = (cos 30 + 3 cos 9) 

= l(2cos40 + 8cos20 + ^^ f 



(33) 
(34) 



If the tune v is far from the resonance lines v = m/n, where m and n are integers, the invariant 
tori of the unperturbed motion are preserved under the presence of the small perturbations by 
virtue of the Kolmogorov-Arnold-Moser (KAM) theorem. If, however, v is close to the above 
mentioned resonance line, under some conditions the KAM invariant tori can be broken. 

Consider first the case where there is only one sextupole located at s = si with f3 x (si). 
Taking the third order resonance, m/3, for example, we keep only the sinusoidal function with 
phase in eq. 3C and the dominant phase independent nonlinear term in eq. 31, and as the 

JT 



result, eqs. 3C and 31 become 



with 



A 



Ji + A sin 3^i 
= ^1 + BJ[ 

{2J 1 (3 x {s 1 )f' 2 (b^L 
P 



b 2 L 



P 



(35) 
(36) 

(37) 
(38) 



where we have dropped the constant phase in eq. [H] and take the maximum value of cos 3 (^I / i), 1. 
It is helpful to transform eqs. [T?] and 38 into the form so-called standard mapping expressed 
as 

J = I + K sin 9 (39) 

9 = 9 + 1 (40) 

with 9 = 3^, / = 3BJ\ and Kq = 3AB. By virtue of the Chirikov criterion ||] it is known 
that when \K \ > 0.97164 p] resonance overlapping occurs which results in particles' stochastic 
motions and diffusion processes. Therefore, 



\Kn\ < 1 



(41) 
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can be taken as a natural criterion for the determination of the dynamic aperture of the machine. 
Putting eqs. ^ and into eq. 41, one gets 

/ 1 7 I T \ 2 

\K \ = 3JM Sl f [—) <1 (42) 
and consequently, one finds maximum J\ corresponding to a m/3 resonance 

^ x (si) 3 \\b 2 \L 
The dynamic aperture of the machine is therefore 



J\ Jmax,sext — oa/„\3\\u\Tl (^3) 



A d yna,sext — \/ 2J m ax,sextPx(s) — /Tin r \a/o ( II Ir ) (^4) 

v y^PxK s \) 1 \\ b i\ L ) 

Eq. 44 gives the dynamic aperture of a sextuple strength determined case. The reader can 
confirm that if we keep sin^i) term instead of sin(3 1 I / i) in eq. 35, one arrives at the same 
expression for A dyna ^ ext as expressed in eq. (Hfl . 

Secondly, we consider the case of single octupole located at s = s 2 with ft x {s 2 ). Taking the 
forth order resonance, m/4, for example, we keep only the sinusoidal function with phase 4*i 
in eq. 3C and the dominant phase independent nonlinear term in eq. [H], and as the result, we 
have eqs. 30 and [H] reduced to 

J~ = J 1 + A sin 4*i (45) 
Wl = + BJl (46) 

with 

a = {JM f? (M) (47) 

B = 2&fe) 2 (^) (48) 

where we have dropped the constant phase in eq. |3l] and take the maximum value of cos 4 (*i), 
1. By using Chirikov criterion, one gets 

Jl < Jmax,oct = 7777^9 ( tAt) ( 49 ) 



W x {s 2 ) 2 \\h\L 

and the corresponding dynamic aperture: 



A dyna , oct = pj max , Ms) = (50) 
v Px\S2) y \03\L 

Thirdly, without repeating, we give directly the dynamic aperture due to a decapole located at 
s = s 3 : 

1 \ 1/6 / p \ 1/3 



A dyna , deca = p(3 x {s) [j^) (51) 

where 64 is the coefficient of the decapole strength. Finally, we give the general expression of 
the dynamic aperture in the horizontal plane [z = 0) of a single 2m (m > 3) pole component: 

1 \2C^2)/ n \ V(m-2) 



where s(2m) is the location of this multipole. Eq. 52 gives us useful scaling laws, such as 

A ( P A 1/(m ~ 2) A A ( 1 \^=T) 

Sidyna,2m OC y\ bm '_^ L J , ana A d yna,2m OC ^ j3 m( s ( 2 m)) ) 

If there is more than one nonlinear component, how can one estimate their collective effect 
? Fortunately, one can distinguish two cases: 
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1) If the components are independent, i.e. there are no special phase and amplitude relations 
between them, the total dynamic aperture can be calculated as: 



Adyna,total — j (53) 

Hi I- Yjj i- J2k i- • • • 

dyna,sext,i dyna,oct,j dyna,deca,k 

If the nonlinear components are dependent, i.e. there are special phase and amplitude 
relations between them (for example, in reality, one use some additional sextupoles to 
cancel the nonlinear effects of the sextupoles used to make chromaticity corrections), there 



is no general formula as eq. 53 to apply. 



In the above discussion we have restricted us to the case where particles are moving in the 



horizontal plane, and the one dimensional dynamic aperture formulae expressed in eqs. 44 



51 and |53|, are the maximum stable horizontal excursion ranges with the vertical displacement 
y = 0. In the following we will show briefly how to estimate the dynamic aperture in 2 dimensions 
when there is coupling between the horizontal and vertical planes. Now we consider the case 
where only one sextupole is located at s = si, and we have the corresponding Hamiltonian 
expressed as follows: 



pi K x {s) 2 p\ K y (s) 2 1 d*B z 
2 ^ 2^2^ 2 y ^ 3\Bp dx 2 

' k=—oo 



H = ^ + ^P* 2 + ? + ^FV + ^^(* 3 - 3xy 2 )L £ 5(s - kL) (54) 



Generally speaking, there exists no universal criterion to determine the start up of stochastic 
motions in 2D. Fortunately, in our specific case, we find out the similarity between the Hamil- 
tonian expressed in eq. ^4] and that of the Henon and Heiles problem which has been much 
studied in literature [11|. The Henon and Heiles problem's Hamiltonian is given by 

Hh^h = \{ x2 +pI+V 2 +P 2 y + V* - (55) 

when Hfjk,H > 1/6 the motion becomes unstable. The intuition we get from this conclusion is 
that there should exist a similar criterion for our problem, i.e. to have stable 2D motion one 



should have H < H max . Note that K x (s) and K y (s) in eq. 54 are equal to unity in the Henon 
and Heiles problem's Hamiltonian. The previous one dimensional result helps us now to find 
H m ax- When y = one has H max oc A 2 dynasextx , since x < A dynaiSext;X . When y ^ 0, the 



crossing terms in eqs. 54 and 55 will play the role of exchanging energy between the two planes, 
and for a given set of x and y the total energy of the coupled system can not exceed H max . If 
we define Ad yn a,sext,y is the dynamic aperture in y-plane, one has: 

Px(S\)A^y na sext x = Py(Sl)A^y na sext y + Px(si)x 2 (56) 

or: 



Adyna,sext,y — J ^ (Ad yna ,sext,x x2 ) ffl) 

where (3 y (si) is the vertical beta function where the sextupole is located and Ad yna ,sext,x is given 
by eq. ^4]. Different from eq. [52], the derivation of eq. 57 is quite intuitive, hinted by the Henon 
and Heiles problem which has been studied numerically instead of analytically in literature. 
From eq. |57] one understands that the difference between Ad yn a,sext,y and Ad yna ,sext,x comes from 
\J I3x{si)/ (5 y (si). If there are many sextupoles in a ring one usually has A dynaiSext>x ^ A dynajSext>y 
since (3 x (si) will not be always larger or smaller than (3 y (si). 
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4 Comparison with simulation results 



To verify the validity of eqs. 44, Mi, R3l and 57, we compare the dynamic apertures of some 



special cases by using these analytical formulae with a computer code called BETA [12]. Taking 
the lattice of Super-ACO as an example, we show the schematic layout of the machine in 
Fig. |. 



BETA— LNS V4.70 /18/09/96/ 7-Oct-99 17:01:44 



1st .00E+00 




saco-full - no sextuple and octupole 



Figure 1: The schematic layout of Super-ACO. 

The horizontal beta function distribution and the working point in the third order tune 
diagram are given in Figs. |2] and ||, respectively. From Fig. [2| one finds that the horizontal 
beta function at the beginning and the end of the one turn mapping is 5.6 m. In the following 
numerical simulations the dynamic apertures correspond to (3 x (s) = f3 x (0) = 5.6 m and it will 
also be used in the analytical formulae. Defining the sextupole, octupole, and decapole strengths 
S = biLj p (1/m 2 ), O = bzL/p (1/m 3 ), D = b^L/p (1/m 4 ), respectively, we make now a rather 
systematic comparison. 



1) 



3) 



4) 



A sextupole is located at s = s± with S(si) = 1 and /3 x (si) = 13.6 m, and its influence 
on the horizontal dynamic aperture is illustrated in Figs. |] and [|. From eq. 44 one gets 
analytically that Ad ynaS ext = 0.0385 m compared with the numerical value of 0.04 m. 

An octupole is located at s = s\ with 0{s\) = 10 and f3 x {si) = 13.6 m, and its influence 
on the horizontal dynamic aperture is illustrated in Figs. || and [?]. From eq. 50 one gets 
analytically that Ad yn a oct = 0.055 m compared with the numerical value of 0.054 m. 



The validity of eq. 51 has been checked also. If a decapole with strength D = 1000 is 
located at s = si, the dynamic aperture in horizontal plane is shown in Figs. | and |. 
From eq. 51 one gets Ady na ,deca = 0.022 m compared with the numerical value of 0.024 m. 



A sextupole of S = 2 and a octupole of O = 62 are located at s = s± and j3 x {s\) = 13.6 
m, and their combined influence on the horizontal dynamic aperture is shown in Figs. [lC| 
From eq. |53] one gets A^ yna tota i = 0.0145 m compared with the numerical value 



and 11 



of 0.016 m. 
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BETA-LHS V4.70 /18/09/96/ 7-Oct-99 17:03:13 



1st .00E+00 



qP^iqAL FlINpTjTONS 




saco-full - no sextuple and octupole 



Figure 2: The horizontal beta function distribution of Super- ACO (@ x (0) = 5.6 m). 



BETA-LHS V4.70 /18/09/96/ 7-Oct-99 17:21:09 

2 NIIZ | | t f : TUNE DfABRAM 



1st .00E+00 




saco-full - no sextuple and octupole 



Figure 3: The tune diagram of the third order of Super-ACO, where the cross indicates the 
working point of the machine. 
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BETA-LHS V4.70 /18/09/96/ 7-Oct-99 14:40:55 

± ! , ! , , , ! , , , ! DYNAMIC APERT^IRE 



1st .00E+00 




0.0 0.02 0.04 0.06 



saco-full - no sextuple and octupole 



Figure 4: The dynamic aperture plot (S(si) = 1 and (3 x (si) = 13.6 m). 



BETA-LHS V4.70 /18/09/96/ 7-Oct-99 14:41:11 

X, XP TRACKING 



1st . 00E+00 



-.06 -.04 



Ex/PI= 1 . 607E 

Ez/PI= 8 . 361E 
Ntour= 100 

Gmax = 1.000E4J05 




02 ft,W «Me-03 Me+00 006 

Z0,ZP0= . 000E+00 . 000E+00 



saco-full - no sextuple and octupole 



Figure 5: The horizontal phase space (S(si) = 1 and f3 x (s\) = 13.6 m). 
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saco-full - no sextuple and octupole 



Figure 6: The dynamic aperture plot (O(si) = 10 and /3 x (si) = 13.6 m). 



BETA— LNS V4.70 /18/09/96/ 7-Oct-99 11:61:11 

X.XP TRACKING 



1st . 00E+00 



Ex/PI= 1.607E-0 

Ez/FI= 8 . 361E^0 

Htour= 100 

Gmax = 1.000E-tj05 



-.04 - 




:+oo 

000E+00 



saco-full - no sextuple and octupole 



Figure 7: The horizontal phase space (O(si) = 10 and P x (si) = 13.6 m). 
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BETA-LHS V4.70 /18/09/96/ 8-Nov-99 15:37:27 



1st .OOE+00 




Ex/PT= 1 . 607E- 
Ez/Pli 8.361E- 
Ntouri 100 
Gmax = 1.000E+05 



saco-full - no sextuple and octupole 



gure 8: The dynamic aperture plot (D(si) = 1000 and (5 x {si) = 13.6 m) 



BETA— LNS V4.70 /18/09/96/ 8-NOV-99 16: II: It 



1st . 00E+00 



-.06 -.04 



Ex/PI= 1 . 607E 

Ez/PI= 8 . 361E 
Ntour= 100 

Gmax = 1.000E-IJ05 




02 &?X90= t.iM-02 Me+00 °- 06 

Z0,ZP0= . 000E+00 . 000E+00 



saco-full - no sextuple and octupole 



ure 9: The horizontal phase space (D(si) = 1000 and (3 x ( s i) = 13.6 m) 
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BETA-LHS V4.70 /18/09/96/ 7-Oct-99 16:39:06 



1st .00E+00 




-.03 -.02 -.01 0.0 0.01 



saco-full - no sextuple and octupole 



Ex/PI= 1.607E- 

Ez/PI= 8.361E- 
Nt.our- 1000 

Gmax = 1.000E-IJ05 



Figure 10: The dynamic aperture plot (S'(si) = 2, O(si) = 62, and f3 x (s\) = 13.6 m 



BETA— LNS V4.70 /18/09/96/ 7-Oct-99 16:41:11 

^ , , , X,XP TpACKIMC 



-Ft: 



1st . 00E+00 



Ex/PI= 1.607E-04 

Ez/PI= 8.361E^09 

Ntour= 100 

Gmax = 1.000E-lj05 



,"*P0= HBOE^oS 2 .000E+00 °- 04 
Z0,ZP0= . 000E+00 . 000E+00 



saco-full - no sextuple and octupole 



Figure 11: The horizontal phase space («S*(si) = 2, 0(s\) = 62, and (3 x (si) = 13.6 m 
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A sextupole of S = 2 and a octupole of O = 62 are located at s = s±, /3 x (si) = 13.6 m, 
and s = S2, 0x( s 2) = 15.18 m, respectively, and their combined influence on the horizontal 
dynamic aperture is shown in Figs. |l2| and |13| . From eq. 53 one gets Ad yn a,totai = 

0.0138 

m compared with the numerical value of 0.0135 m. 

Four sextupoles of S = 2 are located at s = si, 2,3,4 with f3 x {s\) = 13.6 m, (3 x (s2) = 15.18 
m ) Ar( s 3) = 7.8 m, and I3 x {sa) = 6.8 m, respectively, and their combined influence on the 
horizontal dynamic aperture is given in Figs. 14 and 15. From eq. 53 one obtained the 
analytical value A dyna ± oto x = 0.012 m. 

We show how the dynamic apertures depend on the strengths of sextupole and octupole. 



Fig. |16| shows the comparison between the analytical (solid line) and the numerical (dotted) 
results for a sextupole located at s = si in the machine shown in Fig. ^. Obviously, 
Adyna,sext scales with 1/5". Fig. gives the similar comparison for an octupole located 
at s = S2 , and confirms that Ad yna ,oct scales with 1 / \fS. 

Now we change the tune of the machine a little bit (from u x = 1.7 to u x = 1.565), and the 
corresponding horizontal beta function distribution and the third order tune diagram are 
shown in Figs. 18 and 19. It is known that in this case /3 X (0) = 5.1 m. A sextupole of 
S = 2 is located a s = s% with f3 x {s\) = 12.42 m, and its influence on the dynamic aperture 
is shown in Figs. ^ and 21. From eq. 44 one finds Ad yn a,sext = 0.021 m compared with 
the numerical value of 0.02 m. 



Finally, a 2D dynamic aperture is calculated numerically and analytically. If a sextupole 
of S = 2 is located at S2 in the same lattice as in case (1) with (3 X = 15.18 m and (3 y = 4.26 
m, the 2D dynamic aperture is calculated by using BETA and eq. 57 as shown in Fig. 
2| and Fig. 



23 



respectively. The peak analytical dynamic apertures in horizontal and 
vertical planes are 0.0163 m and 0.031 m, compared with the numerical results of 0.017 m 
and 0.034 m, respectively. 
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BETA-LHS V4.70 /18/09/96/ 8-Oct-99 10:31:12 

04 ! , ! , , , ! , , , ! DYNAMIC APERT^IRE 



1st .00E+00 




Ex/PT= 1 
Ez/PI= 8 
NtOur= 1 



6©7E= 
361E- 
000 



saco-full - no sextuple and octupole 



Figure 12: The dynamic aperture plot (S(s\) = 2, 0(s2) = 62, (3 x (s\) = 13.6 m, and (3 x {s-z) 
15.18 m ). 
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Figure 13: The horizontal phase space (S(s\) = 2, 0(s2) = 62, (3 x {s\) = 13.6 m, and f3 x {s2) 
15.18 m). 
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Figure 14: The dynamic aperture plot (£(51,2,3,4) 
Px(s3) = 7.8 m, and /? x (s 4 ) = 6.8 m). 
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BETA— LNS V4.70 /18/09/96/ 10-Oct-99 15:5B:M 



-.06 -.04 



1st . 00E+00 



Ex/PI= 1 . 607Et 

Ez/PI= 8 . 361Et 

Ntour= 100 

Gmax = 1.000E-tj05 



•° 2 &,W .Me-03 Mfc+00 °- 06 

zo,zpo= . 000E+00 . 000E+00 



saco-full - no sextuple and octupole 



Figure 15: The horizontal phase space (£(51,2,3,4) = 2, (3 x (s{) = 13.6 m, (3 x (s2) = 15.18 m, 
Px(s3) = 7.8 m, and P x (s±) = 6.8 m). 
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Figure 16: The dynamic aperture of Super-ACO vs S (S = b^L) p) at s\. 
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Figure 17: The dynamic aperture of Super-ACO vs O [O = b^L/p) at S2- 
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Figure 18: The horizontal beta function distribution of modified Super-ACO ((3 X (0) = 5.1 m). 
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Figure 19: The tune diagram of the third order of modified Super-ACO, where the cross indicates 
the working point of the machine. 
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Figure 20: The dynamic aperture plot (S(s\) = 2 and (3 x (si) = 12.42 m). 
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Figure 21: The horizontal phase space («S*(si) = 2 and (3 x {si) = 12.42 m). 
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Figure 22: The 2D dynamic aperture of Super-ACO with S=2 located at S2 with (5 x {s2) = 15.18 
m and (3 y (s2) = 4.26 m. 
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Figure 23: The analytical estimation of the 2D dynamic aperture of Super-ACO with S=2 
located at S2 with (3 x (s2) = 15.18 m and (3 y {s2) = 4.26 m. 
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To make the comparison much more clear we illustrate the machine parameters in Table 1 
and the comparison results in Table 2. 

The agreement between the analytical and the numerical simulation results is quite good. 
The fact that the analytically estimated 2D dynamic aperture agrees well with that calculated 
by the numerical method shows that the method we have used to treat the two coupled nonlinear 
oscillators is a reasonable heuristic procedure. 



Case 


Multipole strength 


beta function (m) 


1 


S( Sl ) = 2 (1/m 2 ) 


= 13.6 


2 


O(si) = 10 (1/m 3 ) 


P x {s x ) = 13.6 


3 


D(si) = 1000 (1/m 4 ) 


p x ( Sl ) = 13.6 


4 


S( Sl ) = 2 (1/m 2 ), 0{si) = 62 (1/m 3 ) 


P X {8 X ) = 13.6 


5 


S( Sl ) = 2 (1/m 2 ), 0(s 2 ) = 62 (1/m 3 ) 


P x {sx) = 13.6 , 0(s 2 ) = 15.18 


6 


S( Si , 2 ,3a) = 2 (1/m 2 ) 




8 


S( Sl ) = 2 (1/m 2 ) 


= 12.42 , X (O) = 5.1 


9 


S( Sl ) = 2 (1/m 2 ) 


f3 x {s 2 ) = 15.18 



Table 1: summary of parameters 



Case 


Adyna,analy. (m) 


Adyna,numer. (m) 


1 


0.0385 


0.04 


2 


0.055 


0.054 


3 


0.022 


0.024 


4 


0.0145 


0.016 


5 


0.0138 


0.0135 


6 


0.012 


0.0135 


8 


0.021 


0.02 


9 


A x = 0.0163 , A y = 0.031 


A x = 0.017, A y = 0.034 



Table 2: summary of comparison results 



5 Beam-beam interaction limited dynamic apertures 

Since we are interested in the single particle motion under beam-beam forces, the incoherent 
force should be taken into account IfO 



**>-±=^(i— (-£)) m 

where n is the line particle number density, (3 is the particle's velocity in units of the speed of 
light, q is the particle's electric charge, eo is the permitivity in vacuum, r is the transverse offset 
of a particle with respect to the center of the counter-rotating colliding bunch, a is the standard 
deviation of the Gaussian transverse charge density distribution, the positive and the negative 
signs correspond to the colliding bunches with the same or opposite charges, respectively. Now 
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we expand F Ti i n {r) into Taylor series: 

" ^ " ^ r3 + is?'' 5 " ^ + ■ ) <59> 

where we take (3 = 1. To start with, we consider the particle's motion in horizontal plane 
(y = 0) and consider only the delta function nonlinear beam-beam forces coming from one IP. 
The differential equation of motion can be expressed as: 

d 2 x j » nq 2 1 1 o 1 5 

— I +K x {s)x = ± r +^~6 X 

as^ 7reo?n.oc 7 2<r^ 8a 4 48(7° 



, oo 

" 384^ ? + " " " )L S S(s-kL) (60) 

fc=— oo 

where K x (s) describes the linear focusing of the lattice in the horizontal plane, tuqc 2 is the rest 
energy of the particle, 7 is the normalized particle's energy, and L is the circumference of the 
circular collider. The corresponding Hamiltonian is expressed as: 

tt = Px, K *( s ) 2 nq 2 ( 1 2 J_ 4, 1 6 

2 + 2 X + vreomoc^McT 2 ^ 32a 4 288cr 6X 

1 00 
" 3072^ x8 + • • ^ £ (61) 

k=— oo 

where p x = dx/ds. 

To make use of the general dynamic aperture formulae shown in section 3, one needs only 
to find the equivalence relations by comparing two Hamiltonians expressed in eqs. |l| and H 
respectively, and it is found that: 



p 87reomoc 2 7cr 



—L = JVe % A (62) 



= ~ A —z a (63) 

— L = N *\ 8 (64) 

p 3847reomoc 2 7C7 8 

and so on, where we have replaced nL by iV e which is the particle population inside a bunch. 
Till now one can calculate all kinds of dynamic apertures due to nonlinear beam-beam forces. 
For example, one can get the dynamic apertures due to the beam-beam octupole nonlinear force: 



A _ VM^ 

■™-dyna,8,x 



VPJs) ( 87re m c 2 7cr 4 \ 



and 



(65) 



f3y(s 

where s* is the IP position. If we measure dynamic apertures by the beam sizes, one gets: 



Adyna,8,y ~ \ a / *\ (^dyna,8,a; X<2 ) (66) 



/ 2 \ V 2 

A dy na,8,x / 87re m c <ye x \ 
K x , 8 = — "—r^- = 77-^ 67) 
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where e x is the bunch horizontal emittance. When the higher order multipoles effects (2m > 8) 
can be neglected eqs. 55 and ^ give very good approximations to the 2D dynamic apertures 
limited by one beam-beam IP. If there are Njp interaction points in a ring the dynamic apertures 



described in eqs. 35 and 36 will be reduced by a factor of y/Njp (if these Njp interaction points 
can be regarded as independent). Given the dynamic aperture of the ring without the beam- 
beam effect as A XiV , the total dynamic aperture including the beam-beam effect can be estimated 
usually as: 

1 



+ 



oo,x,y 



Taking PEP-II B-Factory design parameters for example []T4| and assuming that the beams 



are round at IP, for the high energy ring, N e = 2.8 x 10 10 , 7 
gets from eq. 67, TZ X ,8 = 3.2, and for the low energy ring, N e 
e x = 49 nm, one finds 1Z X a = 2.7. 



1.76 x 10 4 , and e x = 49 nm, one 
6 x 10 10 , 7 = 6.07 x 10 3 , and 



6 Conclusion 

We have derived the analytical formulae for the dynamic apertures in circular accelerators due to 
single sextupole, single octupole, single decapole (single 2m pole in general), and the combination 
of many independent multipoles. The analytical results have been systematically compared with 
the numerical ones and the agreement is quite satisfactory. These formulae are very useful both 
for the physical insight and in the practical machine design and operation. One application 
of these formulae is to estimate analytically the beam-beam interaction determined dynamic 
aperture in a circular collider. 
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